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Broken  Symmetry  in  Ideal  Magnetohydrodynamic  Turbulence 


John  V.  Shebalin* 

National  Aeronautics  and  Space  Administration 
Langley  Research  Center 
Hampton,  Virginia  23681,  USA 

A  numerical  study  of  the  long-time  evolution  of  a  number  of  cases  of  inviscid,  isotropic,  incompressible, 
three-dimensional  fluid  and  magneto-fluid  turbulence  has  been  completed.  The  results  confirm  that  ideal 
magnetohydrodynamic  turbulence  is  non-ergodic  if  there  is  no  external  magnetic  field  present.  This  is  due 
essentially  to  a  canonical  symmetry  being  broken  in  an  arbitrary  dynamical  representation.  The  broken 
symmetry  manifests  itself  as  a  coherent  structure,  t.e.,  a  non-zero  time-averaged  part  of  the  turbulent 
magnetic  field.  The  coherent  structure  is  observed,  in  one  case,  to  contain  about  eighteen  percent  of  the 
total  energy. 
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1  Introduction 


Several  years  ago,  a  numerical  study  of  isotropic,  inviscid,  incompressible,  two-dimensional  (2D)  turbu¬ 
lence  indicated  that  ideal  magnetohydrodynamic  (MHD)  turbulence  was  non-ergodic  [1].  A  Fourier  spectral 
method  was  used,  and  the  non-ergodicity  manifested  itself  in  the  appearance  of  significant,  non-zero  time- 
averages  for  the  Fourier  coefficients  of  the  turbulent  magnetic  field.  These  results  contrasted  with  the 
predictions  of  canonical  ensemble  theory  [2]  and  indicated  the  presence  of  a  dynamically  broken  symmetry. 

Here,  this  work  is  extended  in  two  ways.  First,  the  dynamics  become  fully  three-dimensional  (3D), 
and  second,  a  more  accurate  (third-order)  time-integration  scheme  is  used  (a  comparison  is  made  with  the 
second-order  scheme  used  previously).  The  results  are  qualitatively  the  same  as  previously  seen,  although  a 
stronger  effect  is  observed.  Also,  the  source  of  the  broken  symmetry  is  more  explicitly  identified. 


2  Basic  Equations 


The  equations  which  describe  ideal,  incompressible  3D  fluid  and  magneto-fluid  dynamics  are 


d,B  =  V  X  (u  X  B) 

=  V  X  (u  X  w-l- j  X  B) 


(1) 

(2) 


Here,  the  fluid  velocity  is  u,  the  vorticity  is  w  =  V  x  u,  the  magnetic  field  is  B  =  b  -I-  Bo  (where  b  is 
the  turbulent  part  and  Bo  is  a  constant,  externally  imposed  part),  and  the  current  is  j  =  V  x  b.  Also, 
Vu  =  VB  =  0. 

If  B  =  0  for  all  time,  then  (1)  is  identically  satisfied  and  (2)  assumes  the  vorticity  form  of  the  incom¬ 
pressible  Euler  equation.  If  Bo  =  0  while  b  ^  0,  then  (1)  and  (2)  are  unchanged,  i.c.,  symmetric,  under  the 
substitution  b  — *  -b.  If,  however,  Bo  ^  0  also,  then  the  equations  are  no  longer  symmetric  under  b  — »  -b. 

Note  that  these  equations  are  invtscid.  It  is  well  known  that  solutions  to  these  equations  are  quite 
different  from  those  for  dissipative  equations,  even  if  the  viscosity  and  resistivity  become  very  small,  so  long 
as  they  are  not  identically  zero.  However,  it  is  possible  to  determine  a  priori  statistical  solutions  in  the 
inviscid  cz«e  (as  will  be  discussed  presently),  while  it  hM  not  proven  possible  to  do  this  in  the  dissipative 
case. 


•  i 
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This  provides  for  a  very  practical  use  of  the  results  presented  herein,  that  is,  as  a  test  of  the  correctness 
of  any  computer  code  which  attempts  to  simulate  incompressible,  homogeneous,  dissipative  turbulence.  In 
any  such  code,  the  viscosity  (and  resistivity)  can  be  set  identically  to  zero;  the  statistical  properties  of  a 
numerical  solution  must  then  be  described  by  the  inviscid  theory  described  here.  There  are,  however,  some 
refinements  which  must  be  incorporated  into  existing  theory;  these  refinements,  the  existence  of  broken 
symmetry  and  non-ergodicity  in  ideal  (3D)  MHD  turbulence,  are  the  subject  of  this  paper. 

3  Canonical  Ensembles 

The  theory  of  canonical  or  ‘absolute  equilibrium’  ensembles,  as  it  applies  to  homogeneous  turbulence,  has 
been  discussed  extensively  before  [1,  3]  (and  the  many  references  therein).  In  brief,  the  physical  variables  u, 
b,  u>,  and  j  are  expanded  in  truncated  Fourier  series,  for  example: 

(3) 

Here,  the  sum  (and  each  similar  sum  appearing  henceforth)  is  over  all  k  such  that  kmin  <  |k|  <  kmai  <  N/2, 
where  N  is  the  number  of  points  in  each  of  the  three  spatial  dimensions.  Since  the  various  fields,  such  as  b(x), 
are  real,  their  coefficients  satisfy  b(k)  =  b*(— k),  where  denotes  complex  conjugation  (here  and  henceforth, 
explicit  time-dependence  will  be  dropped  for  brevity).  Also,  u(k)  =  iib'^k  x  u>(k)  and  j(k)  =  tk  x  b(k), 
and  V  ■  b(x)  =  0  — ►  ik  •  b(k)  =  0,  etc. 

The  independent  real  and  imaginary  parts  of  the  coefficients  u(k)  and  b(k)  cem  be  used  to  label  the  axes 
of  a  multidimensional  phase  space.  The  corresponding  dynamical  system  is  described  by  a  single  point  in 
the  phase  space,  which  moves  about  as  the  system  evolves  in  time.  The  probability  that  the  system  point 
is  in  any  part  of  phase  space  can  be  described  by  a  canonical  distribution  function  which  depends  only  on  a 
small  set  of  conserved  quantities,  the  integral  invariants  of  the  dynamical  system.  Once  the  joint  probability 
distribution  has  been  found,  then  the  equilibrium  energy  spectra  (kinetic  and  magnetic)  are  determined, 
even  though  u(k)  and  b(k)  are  random  variables. 

In  the  case  of  isotropic,  incompressible  3D  Euler  turbulence  (B  =  0),  the  integral  invariants  are  the 
energy  E  and  the  kinetic  helicity  W*.  For  ideal,  isotropic,  incompressible  3D  MHD  turbulence  with  =  0, 
the  integral  invariants  are  the  energy  E,  the  cross  helicity  He,  and  the  magnetic  helicity  Hm\  if  Bo  ^  0,  then 


Hm  >>  no  longer  conaerved,  although  E  and  He  are  still  integral  invariants.  These  three  different  situations 

will  now  be  presented  in  more  detail.  ^ 

S.l  Euler  Turbulence 


Ar 


The  integral  invariants  for  incompressible  3D  Euler  turbulence  [4,  5]  are  ^ 

^=2^*1!  //*  =  2^X;*u(k)  .u.*(k),  (4) 

^Niin 

while  the  distribution  function  is 

/a*  -  ^  * 

D  =  Ctxp(-aE-l3H,)  where  C  =  f]  (5) 

^«isa 

The  product  above  (and  each  product  appearing  henceforth)  is  taken  over  all  k  such  that  kmin  <  |k|  <  kmar. 

but  only  over  independent  values  of  k,  i.e.,  if  k  is  used,  then  — k  is  not.  Using  this  distribution  function,  an  # 

expectation  value  can  be  defined  as 


(Q)  =  /(JofiVuCk),  (6) 

•'  km,. 

with  which  expectation  values  of  the  moments  of  the  Fourier  coefficients  can  be  found: 

(UH(k))  =  {U/(k))  =  0,  (|u«(k)|»)  =  (|u,(k)0  = 

v3/?ir* 

and  (ufi(k)  w/i(k))  =  (u/(k)  w/(k))  =  (7) 

where  the  subscripts  R  and  I  denote  real  and  imaginary  parts,  respectively,  of  the  complex  coefficients,  and 
where  all  the  components  of  u(k)  for  a  given  k  have  equal  expectation  values.  A  straightforward  algebraic 
manipulation  of  these  expectation  values  yields 


a 


Z(^} 

{E){n)-(w*)‘ 


and 

(Q) 


(8) 


where 


•  • 


km.. 


1  ^  I 


^wt  in  ^nsf  n 

Note  that  the  expected  value  of  the  mean  squared  vorticity  (the  enstrophy  D)  is  determined  once  the 
expectation  values  (7)  of  the  ju(k)p  are  known 


3 


Although  the  enstrophy  ft  is  not  a  conserved  quantity,  it  does  have  an  expectation  value,  so  that  the 
parameters  a  and  l3  are  functions  of  only  one  unknown  quantity,  (ft).  The  quantity  Z  is  the  ratio  of  the 
number  of  independent  values  of  k  to  N^,  the  total  number  of  points  in  the  spatial  grid. 


3.2  Ideal  MHD  Turbulence 


The  integral  invariants  of  ideal,  incompressible  3D  MHD  turbulence  with  Bo  =  0  [6]  are 

^  E  “(k)  •  b-(k). 


=  b*(k). 


while  the  distribution  function  is 


D  =  rexp(-o£'- JW, where  C  =  (10) 

(Please  note  that  the  formulas  appearing  in  this  subsection  correct  those  which  appeared  previously  [1].) 

I'slng  this  distribution  function,  expectation  values  (where  (6)  now  includes  integration  over  the  inde¬ 
pendent  components  of  the  b(k))  for  moments  of  the  Fourier  coefficients  can  be  found: 


3  ..3  [  ft  -  V/k') 


(uK(k))  =  (u,(k))  =  o,  (i««(k)0  =  (iu,(k)0  =  5.V3  ;  . 


(b«(k))  =  (b,(k))  =  0,  (|b«(k)|^)  =  (Ibr(k)l^)  =  \ 

*  (o^  —  J^/4)  — 


(u«(k)  b;,(k))  =  (u,(k)  b;(k))  =  ---  '1!'*  2 ■ 

^  -  .v/4)  -  a-*y-7fr^ 


(j«(k)  b-^(k))  =  a,(k)  b;(k))  =  :-.v’  '  ,  (11, 

^  [(ft*  -  ^'/D  -  ft‘7-/t'J 

where  all  the  components  of  ii(k)  for  a  given  k  have  etpial  expectation  values,  and  where  all  the  rompoiKuils 
of  b(k)  for  a  given  k  also  have  eipial  expectation  values  .Note  that  the  energj  spect  ra  iii  the  .VI  HD  case  (II) 
peak  at  low  k  =  |k|.  while  in  the  Fuler  case  (7).  it  peaks  at  high  k- 


4 


A  straightforward  algebraic  manipulation  of  these  expectation  values  yields 


a 


ZR{R+l){E) 


R{Ef  -{R-^  if  {H,f 


^(R+l)m 

R{E) 


7 


iR-^)(E) 

(fl +!)(//„>  ’ 


where 


(12) 


Thus  a,  0,  and  7  are  functions  of  only  one  unknown  quantity,  R,  which  is  the  ratio  of  the  expectation  values 
of  the  magnetic  part  of  the  total  energy  to  kinetic  part  of  the  total  energy. 

The  integral  invariemts  of  ideal,  incompressible  3D  MHD  turbulence  with  B<,  ^  0  are  only  E  and  He,  as 
given  above,  while  Hm  is  no  longer  conserved.  In  this  case,  fl  =  1  and  7  =  0. 


3.3  Distribution  Parameters 

The  parameters  a,  3,  and  7,  which  appear  in  the  canonical  distribution  functions,  must  be  given  specific 
values  in  order  that  definite  predictions  can  be  made.  This  can  be  done  in  a  number  of  ways.  First,  it 
can  be  assumed  that  the  integral  invariants  of  the  various  cases  have  fixed  values  (although  they  actually 
fluctuate  slightly  in  the  time  evolution  of  a  canonical  system).  The  distribution  parameters  then  vary  only 
with  respect  to  a  single  quantity:  (fl)  for  3D  Euler  turbulence  and  R  for  ideal  3D  MHD  turbulence  (except 
for  the  case  where  Bo  ^  0,  where  R  =  1).  The  equilibrium  entropy  [1],  S  -  So  -  logC,  where  C  is  the 
normalizing  coefficient  of  the  distribution  function,  is  thus  also  a  function  of  a  single  quantity.  Since  S  is  a 
minimum  with  respect  to  the  distribution  parameters  [7],  then  this  fact  can  be  used  to  uniquely  determine 
(fl)  or  R,  as  the  case  may  be.  Note  that  this  can  be  done  prior  to  any  numerical  simulation  of  the  time 
evolution  of  an  isotropic,  inviscid  turbulent  system. 

A  second  method  is  to  run  a  numerical  simulation  and  actually  determine  the  time-averaged  values  Q 
of  all  the  quantities  Q  which  are  required,  so  that  Q  can  be  substituted  for  (Q).  In  particular,  fl  — >  (fl)  or 
R  R.  \  third,  and  more  practical  method,  is  to  use  the  time-averaged  values  of  the  necessary  integral 
invariants,  but  determine  (fl)  or  R  through  a  least-squares  procedure.  This  third  method  will  turn  out  to 
be  the  most  useful. 

Once  the  canonical  predictions  of  the  moments  of  u(k)  and  b(k)  are  available,  and  numerical  simulations 
have  produced  corresponding  time-averages,  then  comparisons  can  be  made.  This  is  essentially  a  test  of  the 
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ergodictty  of  the  various  dynamical  systems,  t.e.,  the  equivalence  of  ensemble  and  time  averages.  Numerical 
results,  including  these  time  averages,  are  presented  next,  followed  by  a  discussion  addressing  any  anomalies 
which  are  observed. 

4  Numerical  Results 

As  in  the  previous  work  [1],  a  Fourier  spectral  method  [9]  with  shifted-grid  dealiasing  [10]  was  used  to 
solve  the  dynamical  equations  numerically.  In  the  previous  effort,  however,  the  time-integration  method 
was  a  second-order  Runge-Kutta  scheme  (RK2),  while  the  method  primarily  employed  in  this  work  was  a 
third-order  scheme,  consisting  of  an  Adams- Beishforth  predictor  coupled  with  an  Adams-.Vfoulton  corrector 
{AB."!)  [8].  The  advantages  of  AB3  over  RK‘2  are  higher  accuracy  and  higher  speed,  although  slightly  more 
storage  was  required  in  computer  memory.  In  order  to  better  relate  the  3D  results  obtained  here  to  the  2D 
ones  obtained  previously  [1],  two  3D  runs  which  had  used  AB3  were  partially  re-run  using  RK2.  with  initial 
conditions  remaining  the  same. 

The  numerical  simulations  performed  in  the  course  of  this  work  are  presented  in  Table  1 ,  All  runs  were 
done  on  a  Cray  V.VIP.  using  a  16^  grid  with  k^n^T  =  ^nd  At  =  0.001.  Fach  run  is  designated  by  three 
characters:  the  first  is  either  E  (for  Euler)  or  M  (for  MHD).  th<*  second  is  a  number  representing  a  different 
set  of  initial  conditions,  and  the  last  denotes  the  time-integration  method  used,  A  (for  AB3)  or  R  (for  RK2). 
(The  "initial  conditions’  include  the  values  of  Bo  and  as  well  as  {u>(k).  b(k);  kmm  <  |k|  <  kmar]  at 

t  =  0.)  Each  of  the  integral  invariants  listed  in  Table  1  fluctuated  no  more  than  a  few  parts  per  million 
during  their  respective  simulations. 

In  each  of  the  simulations,  the  initial  spectra  satisfied  (u(k)P  ~  lb(k)p  ~  k^(‘xp(-‘2k'^/kl)  with  ko  =  2 
(although  the  exact  shape  of  the  initial  spectra  was  not  critical;  the  critical  initial  quantities  were  the  values 
of  the  total  energy  and  helicities).  Time-averages  of  the  components  of  u;(k)  and  b(k)  as  well  as  of  their 
squares  (|6TR(lf)P.  R:  real,  etc.)  were  taken,  averaging  every  HO  lime  steps.  (There  are  six  "components' 
of  u;(k)  and  six  of  b(k).  since  the  i—.  y—.  and  c— components  of  lh(«<‘  each  have  a  real  and  an  imaginary 
part.)  The  .second-order  moments  which  correspond  to  the  same  value  of  it  =  |k|  can  be  combined  to  produce 


average  modal  energies 


* 


and  £™(ib)=^^|b(k)p,  (13) 

where  the  sum  ‘®  on^y  over  all  k  with  the  same  value  of  k  =  |k|  and  where  n*  is  the  number  of  wave 
vectors  k  in  such  a  sum.  The  modal  V2tlue8  of  enstrophy  and  mean  squared  current  were  then  Q  =  k^Et(k) 
and  J  =  i^Em(i),  respectively. 

These  modal  energies  can  be  compared  with  those  obtained  from  ensemble-averaging,  once  the  various 
distribution  parameters  are  found,  as  described  previously.  Remember  that  the  distribution  parameters  for 
the  Euler  cases  depended  on  the  average  value  of  the  enstrophy  Q,  while  those  for  the  MHD  cases  depended 
on  the  average  ratio  of  magnetic  to  kinetic  energy  R.  For  example,  in  the  MHD  run  MIA,  a  time-average 
yields  Rt  =  1.50795,  equilibrium  entropy  minimization  gives  R,  =  1.50799,  and  minimizing  the  toot-mean- 
square  (rms)  error  between  time-  and  ensemble-average  modal  energy  spectra  produces  Rrms  =  1.50599. 
Using  the  rms  values  in  all  cases  gives  the  distribution  parameters  presented  in  Table  2. 

To  see  how  well  the  time-  and  ensemble-averaged  modal  energy  spectra  compare,  consider  Figure  1,  where 
representative  spectra  are  shown  (log  =  logic  in  the  Figures).  Figure  la  corresponds  to  rase  .M2A  where 
fcmin  =  1;  the  goodness  of  fit  is  typical  of  all  the  Euler  and  MHD  runs.  Figure  lb  corresponds  to  ca,se  M3A 
where  kmin  =  2;  the  initial  values  of  the  u;(k)  and  b(k)  for  M3A  were  the  same  as  those  of  M2A  except  that 
all  of  the  |k|  =  1  coefficients  were  set  to  zero  and  not  allowed  to  grow  during  the  simulation.  These  graphs 
show  the  average  energy  for  modes  with  the  same  k  =  |k|;  since  there  are  twice  a.s  many  mc.ii.s  with  k  —  2 
as  with  k  =  I,  Figures  la  and  lb  indicate  that  the  modes  with  k  =  kmm  contain  the  same  tola/ energy.  The 
exact  value  of  kmin  is  thus  not  qualitatively  important  in  any  given  simulation. 

It  is  informative  to  consider  the  exact  values  taken  by  the  Fourier  modes  during  their  lime  evolution. 
F-r  ext .  iple.  the  values  of  w(k)  and  b(k)  for  k  =(1.0,01  for  run  MIA  are  shown  in  Figure  2.  where  the  time 
evolution  is  shown  by  plotting  the  real  vs  the  imaginary  parts  of  the  coefficients  for  all  times  betwe>'n  I  =  0 
and  f  =  1000.  The  behavior  is  decidedly  non-ergndic  since  the  components  clearly  do  not  have  zero  mean 
values,  and  thus  do  not  match  the  ensemble  prediction,  ft  is  interesting  to  note  that  although  averagi-  modal 
behavior  (Figure  1)  is  essentially  as  predicted,  the  detailed  evolution  of  the  modal  coefficii-nts  is  not.  ,Mso, 
notice  in  Figure  2  that  we  appear  to  have 

Wj,  ~  tui;  and  65,  ~  »fc;  (H) 
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This  relationship  was  also  apparent  for  all  the  other  coefficients  which  were  examined,  as  long  as  the  coefficient 
mean  values  were  sufficiently  larger  than  the  corresponding  standard  deviation.  The  relationship  (14)  is,  in 
fact,  b  ~  V  X  b,  i.e.,  j  x  b  ~  0.  In  other  words,  the  ideal  MHD  has  become  force-free  [11],  apparently  on  a 
mode-by-mode  basis. 

The  average  coefficients  u>(k)  and  b(k)  can  be  thought  of  as  the  coherent  structure  underlying  these 
simulations  of  ideal,  isotropic  turbulence.  These  coefficients  can  be  used  to  calculate  the  coherent  energy 
present  in  the  turbulence; 

£-:.(t)=^^|b(k)|2,  (15) 

(The  various  helicities  can  also  be  calculated  in  a  similar  manner.)  The  average  coefficients  exist  for  all  cases 
at  t  =  250,  so  that  a  comparison  can  be  made;  the  results  of  this  are  presented  in  Table  3.  Notice  that  the 
coherent  quantities  in  Table  3  are  a  sizeable  fraction  of  the  corresponding  total  quantities  in  Table  1;  for 
example,  the  coherent  energy  in  run  MIA  is  18.8%  of  the  total  energy,  while  the  coherent  magnetic  helicity 
in  the  same  run  is  83%  of  the  total! 

For  those  runs  which  went  beyond  t  =  250,  the  coherent  energies  between  t  —  520  and  t  =  750  are 
shown  in  Figure  3.  It  is  evident  that  the  MHD  simulations  with  Bo  =  0  (MIA  and  M2A)  have  the  largest 
coherent  energies,  and  that  this  coherent  energy  is  primarily  magnetic.  The  MHD  case  with  Bo  =  1  (M4A) 

has  effectively  no  coherent  energy,  while  in  the  Euler  case  (E2A),  the  coherent  energy  is  clearly  decaying 

(although  still  3%  of  the  total  at  I  =  750). 

5  Broken  Symmetry 

In  a  Fourier  representation  (3),  the  basic  equations  (1)  and  (2)  take  the  form 

b(k)  =  ik  X  ^  [u(p)  X  b(q)]  -(- i(k  Bo)  u(k)  (16) 

p+q=k 

u;(k)  =  ik  X  ^  [u(p)  X  w(q) -l-j(p)  X  b(q)]  i(k  •  B„)  j(k)  (17) 

p+q=k 

These  equations  have  the  symmetry  (i.e..  invariance  under)  u»(k)  —  f'^®u>(k).  b(k)  —  f'^®b(k)  which 
merely  reflects  the  isotropy  of  space,  i.e..  place  x  —  x  a  into  (3).  However,  an  average  over  the  displace¬ 
ments  a  merely  reproduces  the  known  fact  that  the  k  =  0  modes  of  w  and  b  are  zero  and  remain  so,  i.e  .  do 
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not  part iripalt*  in  thf’  <lv natnirs  of  th>-  nf'-al  MHIl  l iirt>iiU-n- r 

In  tile  I  a-sf-  wInTr  B  =  n  ih*-  '-(lual  i' ni'  (  l*o  anain  tis..  (mx-  iIi-  -vnnii-'r.  |i  ^  ^ 

!■>  that  all  nif'tiihfr>  of  thf  '•n>»'rnt>l<-  of  (io<>.sit.h  r>-ali/ai io|i»  .f  Vlflh  'url  .iIm  ■  io  '  i  .1 

{u(k).b(k)  km,n  |k|  <  tmar  t  »  powilhk  f^-all/at  Ion  ihi'll  ■»  I*  {ll'k'  f>  k  *.  k  • 

I  hus  what  IS  ralrulatt*d  in  th*-  fiist-riitil*-  avrran*-  l»i  k  1  i-  r»ailv  ti,  k  •  f'  k  '2  I  k  ^ 

h(k)  arui  b.(k)  =  -b(ki  III  whl'h  ruitf  b(  k  I  r  II  aiil'iinati  allv  H  ’Wor  li.  k  >  I  k 
unilffin**!!  if  lhf>  art-  non  /<to  ihfn  th*-  •■nwinlil'  sviimwirv  is  sail  o  t.<  to  4-1.  .1  it  .1  ■  ,  < 

th*"  dvnariiics  ar*-  non-r  rgudtc )  That  ihf’s  ran  h*-  non  /•-r--  an<l  s,il.-.iani  ia!l .  s  l.a-  '  -  .i  in  i  .  . 

th<*  nuivuTiral  rfsulis  prrwnif-d  h»-rfin 


6  Conclusion 

In  this  paper  a  nunienral  study  of  31)  Fdiler  and  nleal  MHI)  iiirhiiPn  *  was  i  r'-s.iio  ■  \  ,11  ■ 

different  cases  were  presented  in  which  the  dynamic  evolution  wa.  allow.-d  !•  pr  ..  ee.i  •  ,  i,-  !■  •  ,■ 

length  of  dimensionless  simulation  time  These  rr-siills  rp-arly  demonsirait  i|i.  pres,  n  .  ,f  •  •  4.-  ,1 

non-ergodicity.  and  coherent  structure  in  ideal  3f)  MHI)  tiirhiilenre 

Realistic  MHD  turbulence  (except  in  superfluids)  is  not  ideal  since  p-sisiive  and  1  iso  ,.|-  pr. .  .  ss.  -  i.s-  ;  ,i. 
energy.  The  results  presented  here  then  address,  rather  than  practical  •as.s  i.nl>  iiumerK.ii  sininUti  ii- 
which  fluid  and  magnetic  dissipation  can  be  set  identically  ecjual  to  zero  lor  iIk-s.-  siiiMilalion-  '<  w-or 
the  ideal  theory  provides  a  very  useful  means  of  testing  computer  codes 

An  interesting,  open  question  is:  Does  the  non-ergodic  behavior  seen  in  the  id’-al  magiiei.  ,h\ dr"d>  nann 
case  also  occur  in  real  geophysical  and  astrophysical  systems  where  dynamo  actoitv  has  )..... n  ,,),s.  r\e.t  I  1 
'' n  answer  to  this  question,  however,  requires  the  introduction  of  dissipation  into  numerical  smnilai  i.  .10 
along  with  a  substantial  increase  in  grid  size.  Although  this  was  beyoinl  the  scop>-  of  the  [.resent  w..rk  it  is 
a  straightforward  extension  of  it.  provided  that  sufficiently  large  coiiq.uters  can  be  utilized 

Addit ionally,  it  should  be  meni ioned  that  there  has  bt^en  some  very  interest  1  hr  wt.rk  in  d>  n.iim.  al  s|. ms 
theory  concerning  broken  symmetry  [12]  A  full  discussion  of  any  connection  to  the  numerical  work  pres,  nted 
here  is  also  beyond  the  scope  of  the  present  work. 

Finally.  I  would  like  to  thank  Dr  M  Y,  llussaini  for  the  op|>ortumty  to  s|iend  a  year  in  residenc..  at 
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ICASE  as  a  Visiting  Scientist,  and  to  thank  Professor  Geoffrey  Lilley  for  his  useful  comments  concerning 
this  paper. 
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Table  1.  Numerical  Simulations 


Run 

Bo 

. 

Total 

integral  Invariants 

■ 

Time 

E 

Hk 

Hr 

//m 

ElA 

0.12 

1 

- 

250 

0.5000 

-.04443 

- 

E2A 

0.12 

1 

- 

750 

0.5000 

3.142 

- 

- 

E2R 

0.25 

1 

- 

250 

0.5000 

3  142 

- 

- 

MIA 

0.35 

1 

0 

1000 

1  000 

- 

0.1326 

0.2129 

MIR 

0.71 

1 

0 

250 

1.000 

- 

0  1326 

0  2129 

M2A 

0.35 

1 

0 

750 

1.000 

0.1326 

0.09973 

M3A 

0.35 

2 

0 

250 

0.9983 

0.1258 

0.09899 

M4A 

0  35 

1 

1 

750 

1.000 

_ : _ 

0  1326 

- 

Table  2.  Distribution  Parameters 


Run 

a 

7 

^rm» 

ElA 

2.624 

0.006825 

- 

17.08 

- 

E2A 

1558 

-2.062 

- 

23  74 

- 

E2R 

15.41 

-2035 

- 

23.79 

- 

MIA 

1.774 

-0.7825 

-1  682 

1  506 

MiR 

1.773 

-0  7824 

-1  681 

- 

1  506 

M2A 

1  547 

-0  7142 

-1.451 

- 

1  206 

M3A 

1  616 

-0  7205 

-2  131 

1  301 

M4A 

1.401 

-0.7075 

1  000 

II 


• • t • • 


I 


Table  3.  Coherency  at  <  =  250 


Run 

m 

m 

mm 

Hk 

ElA 

0.00301 

— 

-.00110 

— 

— 

E2A 

0.0416 

— 

0.305 

— 

— 

E2R 

0.0413 

— 

0.303 

— 

— 

MIA 

0.0880 

0.179 

— 

0.0364 

0.177 

MIR 

0.0919 

0.172 

— 

0.0365 

0.170 

M2A 

0.00562 

0.0656 

— 

0.0158 

0.0635 

M3A 

0.00464 

0.0542 

— 

0.0126 

0.0368 

M4A 

— 

0.000421 

0.0000955 

•  I 
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Figure  1.  Energy  spectra  for  a)  M2A  and  b)  M3A; - theoretical. 


numerical. 


Figure  2.  Components  of  two  vector  modes  of  Ml  A  as  they  evolve  from  t  =  0  to  1000: 
a)  vorticity  and  b)  magnetic  field. 


Figure  3.  Coherent  energies: - kinetic, - magnetic. 
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